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Abstract: 

New  energetic  ionic  liquids  are  investigated  as  potential  high  energy  density  materials. 
Ionic  liquids  are  composed  of  large,  charge-diffuse  cations,  coupled  with  various  (usually 
oxygen  containing)  anions.  In  this  work,  calculations  have  been  performed  on  the 
tetrazolium  cation  with  a  variety  of  substituents.  Density  functional  theory  (DFT)  with 
the  B3LYP  functional,  using  the  6-31  lG(d,p)  basis  set  was  used  to  optimize  geometries. 
Improved  treatment  of  dynamic  electron  correlation  was  obtained  using  second  order 
perturbation  theory  (MP2).  Heats  of  formation  of  the  cation  with  different  substituent 
groups  were  calculated  using  isodesmic  reactions  and  Gaussian-2  calculations  on  the 
reactants.  The  cation  was  paired  with  oxygen  rich  anions  C104 ,  N03",  or  N(N02)2"  and 
those  structures  were  optimized  using  both  DFT  and  MP2.  The  reaction  pathway  for 
proton  transfer  from  the  cation  to  the  anion  was  investigated. 
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Introduction 


There  is  considerable  current  interest  in  ionic  liquids  as  solvents.  No  other  class 
of  solvents  offers  the  versatility  that  ionic  liquids  do:  They  are  typically  thermally  stable, 
have  negligibly  low  vapor  pressure,  high  density,  and  large  liquid  ranges  up  to  400°C. 
Because  of  these  unique  properties,  ionic  liquids  may  be  used  as  electrolytes  for  batteries, 
extraction  media,  and  catalyst  carriers.  The  popularity  of  ionic  liquids  has  also  been 
spurred  on  by  their  classification  as  ‘green,’  due  to  their  negligible  vapor  pressure, 
thereby  decreasing  levels  of  volatile  organic  carbons  in  the  environment  [1], 

l-Butyl-3-methyl-imidazolium  (BMIM)  is  one  of  the  most  widely  used  and  best 
understood  ionic  liquid  cations,  so  when  beginning  to  design  energetic  ionic  liquids,  this 
cation  is  a  natural  place  to  start.  The  unsubstituted  imidazolium  cation  is  shown  in  Figure 
la.  The  energy  content  of  the  ionic  liquid  cation  needs  to  be  raised  for  high  energy 
applications,  and  this  can  be  done  in  several  ways. 

First,  replacing  the  imidazolium  ring  with  a  more  nitrogen  rich  ring  can  increase 
the  energy  content  of  the  ionic  liquid.  Second,  the  hydrogen  or  alkyl  side  chains  can  be 
replaced  with  high  energy  groups,  such  as  -CN,  -NH2,  -N3,  and  -N02.  Finally,  the  anion 
chosen  should  be  oxygen  rich  to  serve  as  an  oxidizer.  Interesting  anions  include  nitrate, 
perchlorate,  and  dinitramide  anions  (Figure  2).  Two  promising  nitrogen-rich  candidates 
are  the  triazolium  (Figure  lb)  and  the  tetrazolium  (Figure  lc)  cations  with  three  and  four 
nitrogens  respectively,  in  the  ring.  These  two  cations  could  present  a  useful  balance 
between  endothermicity  and  thermal  stability.  A  thorough  study  of  the  1 ,2,4-triazolium 
cation  has  recently  been  published  [2].  The  focus  of  the  present  study  is  on  the  electronic 
structure  of  the  tetrazolium  cation. 
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The  energetic  tetrazolium  cations  that  have  been  synthesized  include  1 -amino-4, 5- 
dimethyltetrazolium  [3],  2-amino-4,5-dimethyltetrazolium  [3],  2,4,5-trimethyltetrazolium 
[4],  and  l,5-diamino-4-H-tetrazolium  [5].  Melting  points  for  these  cations  combined  with 
iodide,  nitrate,  and  perchlorate  anions  range  from  -59  °C  to  156  °C  and  their  liquid  range 
can  be  up  to  229  °C.  A  summary  of  energetic  tetrazolium  cations  that  have  recently 
appeared  in  the  literature  and  their  melting  points  are  given  in  Table  1  [3-6], 

Synthesis  of  triazolium  cations  is  easier  than  that  of  tetrazolium  cations,  hence, 
there  have  been  many  more  studies  of  energetic  triazolium  than  tetrazolium  compounds. 
Triazolium  cations  have  been  successfully  substituted  with  azido,  diazido  [4],  and  amino 
groups  [3-6],  The  successful  synthesis  of  the  triazolium  cation  is  not  so  surprising  due  to 
its  similarity  to  the  popular  imidazolium  cation,  but  what  is  surprising  is  the  use  of  the 
triazole  ring  as  an  ionic  liquid  anion.  Ionic  liquids  typically  contain  a  large  asymmetric 
organic  cation,  which  causes  the  ions  to  be  poorly  coordinated,  but  as  shown  by  Katritzky 
et  al.,  [7]  the  anion  can  be  large  as  well.  An  ionic  liquid,  l-butyl-3-methylimidazolium 
3,5-dinitro-l,2,4-triazolate,  containing  a  high  energy  planar  anion,  has  been  synthesized 
and  has  a  low  melting  point  of  35  °C.  Delocalization  of  charge  on  the  anion  ring  is 
caused  by  the  nitro  substituent  groups,  which  also  serve  to  increase  the  energy  of  the 
compound. 

Interactions  in  ionic  liquids  are  more  complicated  than  those  in  simpler  liquids, 
making  them  more  difficult  to  understand  on  a  molecular  level.  Theory  can  provide  an 
excellent  tool  for  understanding  the  structure  and  dynamics  of  ionic  liquids.  There  have 
been  several  molecular  dynamics  simulations  of  ionic  liquids,  in  order  to  model  bulk 
properties,  such  as  melting  points,  diffusion,  and  viscosity  [8], [9].  Radial  distribution 
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functions  and  densities  have  also  been  calculated  [8], [10].  Predicted  equilibrium 
structures  of  l-alkyl-3-methylimidazolium  ionic  liquids  with  various  anions  and  alkyl 
side  chains  have  been  found  to  compare  well  to  experiment  [11].  Dynamics  simulations 
have  revealed  that  l-ethyl-3-methylimidazolium  nitrate  has  diffusion  properties  similar  to 
those  of  a  supercooled  liquid  [12], 

The  solvent  properties  of  ionic  liquids  and  their  propensity  to  be  synthesized 
depend  greatly  on  their  acid-base  properties,  so  the  ability  to  predict  these  properties 
would  be  of  great  use.  The  more  acidic  the  cation,  the  more  difficult  it  will  be  to 
protonate  the  neutral  ring  to  form  an  ionic  liquid.  Of  the  three  neutral  rings,  imidazole, 
triazole,  and  tetrazole,  the  weakest  base  is  tetrazole.  Tetrazolium  is  therefore,  the  most 
difficult  to  synthesize.  Acidity  constants  of  triazole  and  tetrazole  have  been  calculated 
using  several  semi-empirical  methods  [13], 

In  order  to  have  low  melting  ionic  liquids,  the  charge  of  at  least  one  of  the  ions 
must  be  delocalized.  MCSCF  analysis  of  the  triazolium  cation  shows  that  the  electrons 
on  the  cation  are  shared  between  two  resonance  structures  [2],  This  study  also 
investigated  the  effects  of  more  energetic  and  less  energetic  substituent  groups.  For 
example,  the  nitrile  group  was  proposed  as  a  better  substituent  for  high  energy 
applications  than  an  azide  group.  The  structures  of  triazolium  dinitramide  systems  were 
investigated  using  dimer  pairs.  A  wide  variety  of  geometries  were  found  and  the 
presence  of  small  barriers  for  proton  transfer  from  the  cation  to  the  anion  show  that 
deprotonation  may  be  an  important  mechanism  in  decomposition  of  triazolium-based 
ionic  liquids. 
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Computational  Methods 

Initial  structures  were  obtained  by  performing  gas  phase  density  functional  theory 
(DFT)  calculations  on  the  isolated  ions  using  the  Becke  three  parameter  Lee- Yang-Parr 
hybrid  functional  (B3LYP)  [14,  15].  The  basis  sets  used  were  6-31G(d,p)  [16-18],  and  6- 
31  lG(d,p)  [18,  19].  The  DFT  geometries  and  energies  are  compared  to  second  order 
Moller-Plesset  perturbation  theory  (MP2)  calculations  [20,  21],  A  MCSCF  population 
analysis  [22]  using  Edmiston-Ruedenberg  type  localized  molecular  orbitals  (LMO)  [23] 
was  carried  out  to  investigate  the  amount  of  electron  delocalization  in  the  ring.  The 
orbitals  included  in  the  MCSCF  active  space  are  both  jt  orbitals,  their  corresponding 
antibonds,  and  the  lone  pair. 

Calculations  on  ion  pairs  provide  information  about  the  fundamental  interactions 
between  the  cation  and  the  anion.  The  gas  phase  ion  interactions  can  provide  insight  into 
the  bulk  liquid  structure.  Dimer  pairs  were  optimized  using  DFT  and  MP2  methods  and 
the  6-31+G(d)16  [24,  25]  basis  set.  Hessians  (matrices  of  energy  second  derivatives)  are 
used  to  determine  if  stationary  points  are  minima  or  transition  states.  At  the  final  MP2/6- 
31  lG(d,p)  geometries,  improved  relative  energies  were  obtained  for  some  of  the 
tetrazolium  cation  isomers,  using  singles  and  doubles  coupled  cluster  theory  with 
perturbative  triples  (CCSD(T))  [26,  27]  with  the  6-31  lG(d,p)  and  cc-pVTZ  [28]  basis 
sets. 

Intrinsic  reaction  coordinate  (IRC)  calculations  were  used  to  connect  transition 
states  with  reactants  and  products  [29].  The  step  size  used  for  the  IRC  calculations  was 
0.05  (amu)1/2  bohr.  All  calculations  were  done  with  GAMESS  [30,  31],  and  all  molecules 
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were  visualized  with  MacMolPlt  [32], 


Results  and  Discussion 

Relative  Energies  of  Tetrazolium  Cations 

The  tetrazolium  cation  has  four  parent  isomers  (all  R=H)  labeled  I,  II,  III,  and  IV 
(Figure  3).  Each  isomer  has  two  possible  resonance  contributors.  The  relative  energies 
of  I,  II,  III,  and  IV  are  shown  in  Table  2.  Isomer  I  is  predicted  by  MP2  to  be  the  lowest  in 
energy,  with  isomer  II  only  1.9  kcal/mol  higher  in  energy.  The  MP2  relative  energy 
ordering  is  I<II<IV<IIL  In  general,  DFT  and  MP2  are  in  good  agreement,  although  DFT 
slightly  reverses  the  order  of  isomers  I  and  II. 

The  relative  energies  of  isomers  I  and  II  were  also  calculated  using  the  CCSD(T) 
method  at  the  MP2/6-31  lG(d,p)  geometries.  Using  the  6-31  lG(d,p)  and  cc-pVTZ  basis 
sets,  isomer  II  is  predicted  to  be  lower  in  energy  than  I  by  0.2  and  0.5  kcal/mol 
respectively.  The  open  chain  form  of  the  cation,  azido  formidinium  (V),  shown  in  Figure 
4,  was  also  compared  to  the  four  parent  isomers.  At  the  MP2/6-31 1++G(d,p)  geometry 
this  isomer  is  0.7  kcal/mol  lower  in  energy  than  isomer  I  (Table  3).  According  to 
CCSD(T)/cc-PVTZ//MP2/6-31+G(d),  isomer  V  is  1.8  kcal/mol  higher  in  energy  than 
isomer  I,  and  1.3  kcal/mol  higher  in  energy  than  isomer  II. 

The  large  MP2/6-31  lG(d,p)  barriers  for  proton  transfer  reactions  between  isomer 
I  and  isomers  II,  III,  and  IV  (Figure  5)  imply  that  movement  of  a  proton  on  the  ring  is 
not  likely. 

Calculated  MP2/6-31  lG(d,p)  _  and  it  bond  orders  [33]  (Figure  6)  suggest  that  the 
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two  double  bonds  in  the  ring  are  delocalized.  The  bond  lengths  between  atoms  in  the 
rings  do  not  exhibit  any  significant  changes  among  the  four  parent  isomers  (Table  4).  The 
bond  lengths  between  atoms  in  the  cation  rings  are  consistent  with  a  delocalized  ring 
system  (cf.,  1.32  A  for  a  CN  double  bond  and  1.34  A  for  a  NN  double  bond). 

Bond  lengths  can  suggest  which  resonance  structures  are  favored  for  isomers  I 
and  III.  (The  two  resonance  structures  for  isomers  II  and  IV  are  equivalent.)  The  nearly 
equal  N1-N2  and  N2-N3  bond  lengths  in  isomer  I  (1.31  and  1.29  A,  respectively)  suggest 
that  resonance  structures  I  and  I'  make  similar  contributions  to  the  electronic  structure  of 
the  ring,  with  perhaps  a  slight  preference  for  I  over  I'.  The  two  resonance  structures  for 
isomer  III  are  likewise  equally  important  as  indicated  by  the  equal  N2-N3  and  N3-N4 
bond  lengths  and  the  nearly  equal  lengths  of  the  N1-C5  and  N4-C5  bonds. 

MCSCF  jt  LMO  populations  and  jt  bond  orders  (Table  5)  can  provide  a  more 
sophisticated  analysis  of  competing  resonance  structures  [2].  The  diagonal  density  matrix 
elements  give  the  electron  occupancy  of  the  localized  molecular  orbitals.  The  off- 
diagonal  elements  give  the  bond  orders.  A  positive  bond  order  indicates  a  bonding 
interaction  and  a  negative  bond  order  indicates  an  antibonding  interaction  between  the 
LMOs. 

For  isomer  I,  the  jt  orbital  populations  on  N1  and  N3  are  1.42  and  1.37 
respectively.  In  a  resonance  structure  that  was  purely  I,  N1  would  have  a  jt  population  of 
2.00  and  all  other  atoms  in  the  ring  would  have  jt  populations  of  1.00  (see  Figure  3).  If 
the  structure  was  purely  I’,  then  N3  would  have  a  jt  population  of  2.00  and  all  other 
atoms  in  the  ring  would  have  jt  populations  of  1.00  (cf.,  Figure  3).  In  reality  both  N1  and 
N3  have  jt  populations  of  less  than  1.5  electrons.  This  indicates  that  the  ring  is  a  hybrid 
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between  I  and  I’  and  is  delocalized.  The  difference  in  it  orbital  populations  indicate  that 
I  is  slightly  favored  over  I’.  In  both  I  and  I’  C5-N4  is  a  double  bond.  This  bond  has  a  it 
bond  order  of  0.73,  the  largest  of  the  five  bonds.  The  N2-N3  it  bond  order,  a  double  bond 
in  I,  is  0.65.  The  N1-N2  it  bond  order,  is  a  double  bond  in  I’  is  0.57.  Neither  of  these 
bond  orders  is  close  to  1 .00  indicating  once  again  that  the  ring  is  delocalized. 

A  similar  analysis  can  be  done  for  isomers  II-IV  showing  that  all  the  it  electrons 
are  delocalized  on  the  cations.  According  to  the  it  orbital  populations  in  isomer  III,  the 
lone  pair  is  distributed  equally  on  N1  and  N2.  Similar  comments  apply  to  isomer  II(IV), 
in  which  the  lone  pair  is  distributed  equally  on  N1  and  N4  (N2  and  N3). 

The  relative  energies  of  the  parent  isomers  with  a  single  substituent  on  a  nitrogen 
or  a  carbon  are  shown  in  Table  6.  The  isomers  with  a  substituent  R’  (on  a  nitrogen)  are 
labeled  I-A  though  IV-A  (Figure  7a),  the  isomers  with  a  substituent  R”  (on  the  carbon) 
are  labeled  I-B  though  IV-B  (Figure  7b),  and  the  isomers  with  a  substituent  R”’  (on  a 
nitrogen)  are  labeled  I-C  through  IV-C  (Figure  7c).  For  all  substituents,  the  relative 
energies  remain  the  same  as  those  of  the  parent  isomers  with  no  substitution. 

A  -N02  substituent  does  not  bind  well  to  a  nitrogen  of  the  tetrazolium  ring:  the 
NN  distance  is  2.4  A,  compared  with  a  normal  single  NN  bond  distance  of  1.430  A  as, 
for  example,  in  N2H4  [34],  Because  of  the  lack  of  a  strong  NN  bond,  the  nitro  group 
easily  changes  positions  on  the  ring.  For  example,  III-A(N02)  converts  to  isomer  I- 
A(N02),  with  no  barrier.  A  cation  with  a  -N02  substitutent  bound  to  a  ring  nitrogen,  is 
most  likely  an  ion-dipole  complex,  in  which  the  cation  becomes  a  neutral  tetrazole  ring 
and  the  nitro  group  becomes  an  incipient  cation.  This  is  easily  verified  by  examining  the 
charge  on  the  substituent  -NOz  group.  The  average  of  the  Mulliken  charges  [35]  on  the 
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-N02  group  for  all  isomers  with  -N02  on  N  is  +0.83.  The  average  MP2/6-31  lG(d,p) 
binding  energy  of  N02+  to  a  tetrazole  ring  is  22.8  kcal/mol. 

A  comparison  of  relative  energies  (averaged  for  all  four  parent  isomers)  for 
substitution  on  a  C  vs.  N  is  shown  in  Table  7.  Substitution  at  C  is  almost  always 
energetically  favored  over  substitution  at  N.  The  only  exception  is  for  the  -N02 
substituent.  The  bond  lengths  within  the  ring  do  not  change  significantly  from  the  parent 
isomer  when  a  single  substituent  is  present. 

Heats  of  Formation 

The  energy  content  of  a  heterocycle  can  be  raised  dramatically  by  increasing  the 
number  of  nitrogen  atoms  in  the  ring.  For  example,  the  experimental  heats  of  formation 
for  imidazole,  1,2,4-triazole,  and  tetrazole  are  14.0,  26.1,  and  56.7  kcal/mol  respectively 
[33], 

In  this  work,  heats  of  formation  for  the  cation  rings  were  calculated  using 
isodesmic  reactions  and  the  Gaussian-2  (G2)  method  [36,  37],  a  multilevel  method 
designed  to  obtain  accurate  thermochemistry.  The  G2  calculations  were  performed  using 
Gaussian  94  [38],  Isodesmic  reactions,  in  which  the  number  of  formal  bond  types  is 
conserved,  minimize  the  change  in  correlation  energy  [39],  thereby  reducing  the  error  in 
a  computed  heat  of  formation. 

The  isodesmic  reactions  of  each  resonance  structure  are  shown  in  Figure  8.  For  I 
and  III  each  resonance  structure  has  a  different  isodesmic  reaction.  The  heats  of 
formation  for  I  and  I’  (III  and  III’)  are  similar,  and  their  average  will  be  taken  as  a  heat  of 
formation  for  the  composite  I/F  (III/IIF)  [2],  The  largest  difference  in  heats  of  formation 
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between  resonance  structures  is  13.0  kcal/mol  for  the  -N02  substituent;  the  other 
differences  are  all  less  than  3  kcal/mol. 

The  heats  of  formation  for  the  parent  isomers  are  shown  in  column  2  of  Table  8. 
These  compare  well  to  heats  of  formation  that  were  calculated  previously  using 
Gaussian-3  [13].  The  changes  in  heats  of  formation  for  substituted  cations  relative  to  the 
hydrogen  substituted  cations  are  shown  in  Table  9.  For  the  most  part,  substitution 
destabilizes  the  ring,  although,  substitution  by  a  fluorine  or  an  amine  on  a  carbon 
stabilizes  the  ring  by  up  to  40  kcal/mol.  However,  an  azide  group  or  a  nitrile  group  can 
increase  the  heat  of  formation  of  the  ring  by  as  much  as  1 12  kcal/mol.  Substitution  with 
-F,  -NH2,  and  -N02  generally  increases  the  heats  of  formation  by  about  20  kcal/mol. 

Proton  Transfer  Reactions 

One  possible  reaction  pathway  for  the  simple  ion  pairs  is  proton  transfer  from  the 
cation  to  a  partner  anion  to  form  a  neutral  pair.  The  heats  for  the  reaction  CN4H3+  +  Y 
CN4H2  +  HY  (where  Y=N03,  C104,  or  N(N02)2)  are  shown  in  Table  10.  All  of  the  energy 
differences  (AE)  are  exothermic,  by  46-144  kcal/mol,  and  the  AEs  for  proton  transfer  to 
each  of  the  three  anions  are  similar.  The  heat  of  reaction  for  removal  of  H”\  shown  in 
column  three  of  Table  10,  is  more  exothermic  than  the  removal  of  H’  (see  Figure  3) 
because  the  2,5-disubstituted  neutral  ring  is  3.7  kcal/mol  lower  in  energy  than  the  1,5- 
disubstitued  neutral  ring.  The  1,3-disubstituted  ring  is  considerably  higher  in  energy. 
There  are  three  possible  positions  for  proton  transfer  to  the  dinitramide  anion,  but 
dinitramic  acid  is  most  stable  with  the  proton  on  the  center  nitrogen,  so  this  is  the  acid 
used  when  calculating  the  AEs. 
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Ion  Pair  Interactions 


Ionic  liquid  dimer  pairs  were  optimized  with  DFT  and  MP2  methods  using  the  6- 
31+G(d)  basis  set.  ZPE  corrections  at  the  B3LYP/6-31+G(d)  level  of  theory  were  used  in 
all  cases.  The  structures  of  the  tetrazolium  dimers  with  perchlorate  and  nitrate  partners 
are  shown  in  Figures  9  and  10,  respectively.  For  the  cation  with  perchlorate  or  nitrate,  at 
both  levels  of  theory  a  proton  transfers  from  the  cation  to  the  anion  during  geometry 
optimization  of  the  ion  pair.  Hydrogen  bonding  and  other  attractive  electrostatic 
interactions  are  found  to  be  important  in  the  dimer  structures.  Optimizations  give 
geometries  in  which  the  number  of  attractive  electrostatic  interactions  is  maximized.  For 
pairs  containing  HN03  or  HC104,  there  is  always  one  linear  O-H — N  hydrogen  bond. 
Hydrogen  bond  lengths  vary  from  1.795  A  to  1.846  A  for  HN03  and  1.771  A  to  1.822  A 
for  HC104.  The  majority  of  structures  also  contain  a  non-linear  N-H — O  hydrogen  bond 
or  a  C-H — O  attractive  electrostatic  interaction. 

Tetrazolium  dinitramide  (Figure  11)  pairs  have  a  much  wider  variety  of  structures 
than  tetrazolium  nitrate  or  tetrazolium  perchlorate  pairs.  This  is  because  the  dinitramide 
anion  has  the  ability  to  bind  a  proton  at  three  different  positions  (Figure  12).  For 
HN(N02)2  pairs,  hydrogen  bond  lengths  vary  from  1.721  A  to  1.942  A.  Once  again 
hydrogen  bonding  and  other  electrostatic  interactions  are  important  for  structures 
containing  HN(N02)2.  However,  six  out  of  19  structures  do  not  contain  a  linear  hydrogen 
bond  at  all.  These  structures  usually  contain  weaker,  nonlinear,  N-H — N,  O-H— N,  or  N- 
H — O  attractive  electrostatic  interactions. 

The  relative  energies  of  the  pairs  with  HN03  or  HC104  show  that  the  lowest 


This  nape  is  Distribution  A:  annroved  for  nublic  release:  distribution  unlimited. 


11 


energy  pairs  contain  the  1,3-H  substituted  tetrazole  isomer.  For  example,  structure  1  in 
Figure  11  is  2.8  kcal/mol  lower  in  energy  than  structure  6.  This  difference  in  energy  can 
be  attributed  to  the  difference  in  energy  of  the  1,3-H  substituted  tetrazole  ring  (number  1) 
and  the  1,2-H  substituted  tetrazole  ring  (number  6).  The  lowest  energy  structures  with 
HN(N02)2,  contain  the  1,3-H  substituted  tetrazole  isomer  and  have  the  proton  transferred 
to  the  central  nitrogen  of  N(N02)2".  Also  shown  in  Figures  9,  10,  and  1 1  are  the  changes 
in  energy  for  the  reaction  in  which  the  infinitely  separated  anion  and  cation  come 
together  to  form  either  an  ion  pair  or  a  neutral  pair  by  proton  transfer  from  the  cation  to 
the  anion: 

[Tetrazolium]+  +  [Anion]"  ->  [Tetrazole]  [H- Anion]  or 
[Tetrazolium]+ +  [Anion]'  ->  [Tetrazolium]  [Anion] 

(Anion  =  NOs,  C104,  or  N(N02)2). 

For  pairs  in  which  a  proton  transfers,  the  neutral  pair  is  lower  in  energy  than  the 
separated  neutral  products,  due  to  the  stabilization  from  electrostatic  interactions  in  the 
neutral  dimer  products.  The  average  dimer  stabilization  is  7.3,  9.6,  and  8.4  kcal/mol  for 
[HN03],  [HC104  ],  and  [HN(N02)2]  respectively. 

At  the  DFT  level  of  theory  there  are  some  stable  HN(N02)2"-cation  pair 
geometries;  however  when  those  DFT  geometries  are  optimized  with  MP2,  the  acidic 
proton  usually  transfers  to  the  anion.  Only  one  stable  anion-cation  pair  was  found  at  the 
MP2  level  of  theory.  This  ion  pair  is  structure  13  in  figure  11.  The  ion  pair  is  higher  in 
energy  by  18.9  kcal/mol  than  the  lowest  energy  neutral  pair  at  MP2/6-31+G(d).  In  this 
structure  the  N-H  bond  length  is  stretched  to  1.076  A  as  compared  to  1.019  A  for  the 
non-acidic  hydrogen.  The  hydrogen  bond  length  is  also  much  shorter  at  1.598  A.  The 
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length  of  the  hydrogen  bonds  show  that  the  hydrogen  bond  in  the  ionic  structure  is  much 
stronger  than  the  hydrogen  bond  in  the  various  neutral  structures.  The  MP2  barrier  for 
proton  transfer  from  the  cation  to  the  anion  of  structure  13  of  Figure  1 1  is  0.4  kcal/mol 
(Figure  13),  however  when  ZPE  corrections  are  accounted  for,  the  barrier  disappears. 

For  all  dimer  pairs,  proton  transfer  from  a  cation  to  the  nitrogen  of  the  dinitramide  anion 
is  favored  energetically  over  proton  transfer  to  one  of  the  oxygens  of  the  dinitramide 
anion. 

Gas  phase  calculations  on  an  ionic  pair  can,  of  course,  only  give  an  approximation 
to  the  true  interactions  in  a  crystalized  ionic  liquid.  Proton  transfer  from  the  cation  to  the 
anion  should  be  less  likely  in  the  crystal  or  in  bulk  liquid.  For  example,  in  studies  of 
ammonium  salts,  proton  transfer  occurs  in  the  isolated  gas  phase  ion  pair,  but  stable  ion 
pairs  are  found  when  two  anions  and  two  cations  are  present  [40],  MP2/6-31++G(d,p) 
calculations  on  1 ,2,4-triazolium  dinitramide  typically  do  not  give  stable  ion  pairs  [2], 

The  most  recent  MP2  calculations  done  on  three  ion  pairs  of  1,2,4-triazolium  dinitramide, 
show  that  the  6  ion  cluster  is  only  1.0  kcal/mol  higher  in  energy  than  the  6  neutral  cluster 
[41].  Whether  the  ionic  form  or  the  neutral  form  is  lower  in  energy  for  two  ion  pairs  is 
still  under  investigation.  Future  studies  on  larger  clusters  of  ionic  pairs  are  needed  to 
more  fully  understand  the  interactions  in  the  bulk  liquid. 

Conclusion 

In  this  work,  tetrazolium  cations  I  and  II  were  established  to  be  the  lowest  energy 
isomers,  for  the  parent  isomers  as  well  as  for  all  substituted  isomers.  The  DFT  and  MP2 
energies  are  generally  in  good  agreement.  The  relative  energy  of  the  open  chain  form  of 
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the  cation  is  predicted  to  be  only  slightly  higher  in  energy  than  isomers  I  and  II.  A 
MCSCF  jt  orbital  analysis  indicates  that  the  electrons  in  the  cation  ring  are  delocalized. 
Calculated  heats  of  formation  show  that  the  tetrazolium  cation  ring  has  the  potential  to 
release  large  amounts  of  energy  during  decomposition  and  thus  has  excellent  potential  as 
a  high  energy  fuel.  This  is  especially  true  when  the  ring  is  substituted  with  -N3  or  -CN. 
When  a  cation  is  paired  with  oxygen  rich  anions,  a  single  gas  phase  ion  pair  was  not 
generally  found  to  be  stable.  A  proton  transfers  without  barrier  from  the  cation  to  the 
anion  to  form  a  neutral  pair. 
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Table  1:  Energetic  tetrazolium  cations,  their  melting  points  (Tm),  thermal  decomposition 
temperatures  (Td),  and  glass  transition  temperature  (Tg). 


1  -  Amino-4, 5-dimethyltetrazolium  Iodide3 

Tm  121 °C 

1  -  Amino-4, 5-dimethyltetrazolium  Nitrate3 

Tg  -59  °C 

Td  170  °C 

1  -  Amino-4, 5-dimethyltetrazolium  Perchlorate3 

Tm  51  °C 

Td  182  °C 

2- Amino-4, 5-dimethyltetrazolium  Iodide3 

Tm  124  °C 

2- Amino-4, 5-dimethyltetrazolium  Nitrate3 

Tm  94  °C 

Td  173  °C 

2- Amino-4, 5-dimethyltetrazolium  Perchlorate3 

Tm  140  °C 

Td  238  °C 

2,4,5-Trimethyltetrazolium  Iodideb 

Tm  156  °C 

2,4,5-Trimethyltetrazolium  Nitrateb 

Tm  94  °C 

Td  193  °C 

2,4,5-Trimethyltetrazolium  Perchlorateb 

Tm  133  °C 

Td  315  °C 

1 ,5-Diamino-4-H-Tetrazolium  Perchlorate0 

Tm  125-130  °C 

3  from  Xue,  Twamley,  and  Schreeve3, b  from  Xue,  Gao,  Twamley,  and  Schreeve4, e  from  Drake, 
Hawkins,  Boatz,  Hall,  Vij 5 


Table  2:  Relative  energies  (kcal/mol)  of  Isomers  I,  II,  III,  and  IV. 


Cation 

B3LYP/6-3 1  G(d,p) 

B3LYP/6-31  lG(d,p) 

MP2/6-31G(d,p) 

MP2/6-311G(d,p) 

I 

0 

0 

0 

0 

II 

-0.6 

-0.4 

1.4 

1.9 

III 

16 

16.2 

19.1 

19.2 

IV 

14.9 

14.7 

15.2 

14.7 

Table  3:  Relative  energies  (kcal/mol)  of  the  closed  chain  isomers  I,  II,  and  the  open  chain  isomer,  V, 
using  CCSD(T). 


MP2/6-31 1++G(d,p) 

CCSD(T)/cc-PVTZ// 
MP2/6-31 1++G(d,p) 

CCSD(T)/6-31  lG(d,p)// 
MP2/6-3 1 1++G(d,p) 

I 

0.0 

0.0 

0.0 

II 

1.9 

-0.5 

-0.2 

V 

-0.7 

1.8 

1.6 

This  nape  is  Distribution  A:  annroved  for  nublic  release:  distribution  unlimited. 


17 


Table  4:  Bond  lengths  (Angstroms)  in  the  ring  for  isomers  I  through  IV. 


C5-N1 

N1-N2 

N2-N3 

N3-N4 

N4-C5 

I  1.36 

1.31 

1.29 

1.33 

1.31 

II  1.33 

1.35 

1.30 

1.35 

1.33 

III  1.35 

1.32 

1.32 

1.32 

1.34 

IV  1.35 

1.31 

1.32 

1.31 

1.35 

Table  5:  MCSCF  n  orbital  populations  and  bond  orders. 

jt  orbital  populations 

N1 

N2 

N3 

N4 

C5 

Total  bonding 

I 

1.42 

1.10 

1.37 

1.12 

0.98 

1.93 

II 

1.44 

1.09 

1.09 

1.44 

0.94 

1.92 

III 

1.46 

1.46 

1.04 

1.10 

0.95 

1.91 

IV 

1.07 

1.43 

1.43 

1.07 

0.99 

1.94 

Adjacent  jt  bond  orders 


C5-N1 

N1-N2 

N2-N3 

N3-N4 

N4-N5 

I 

0.55 

0.57 

0.65 

0.48 

0.73 

II 

0.65 

0.44 

0.77 

0.44 

0.65 

III 

0.61 

0.44 

0.60 

0.62 

0.64 

IV 

0.64 

0.56 

0.50 

0.56 

0.64 

Next  neighbor  jt  bond  orders  (antibonding) 

C5-N2 

N1-N3 

N2-N4 

N3-C5 

N4-N1 

I 

-0.22 

-0.30 

-0.27 

-0.11 

-0.15 

II 

-0.18 

-0.19 

-0.19 

-0.18 

-0.29 

III 

-0.30 

-0.33 

-0.22 

-0.15 

IV 

-0.14 

-0.34 

-0.34 

-0.14 
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Table  6:  MP2/6-311G(d,p)  relative  energies  (kcal/mol)  isomers  with  a  single  substitution  on  a 
nitrogen  or  a  carbon.  (See  Figure  7  for  notation.) 


R'= 

I-A 

II-A 

III-A 

IV- A 

H 

0.0 

1.9 

19.1 

14.7 

F 

0.0 

3.4 

18.3 

12.6 

CN 

0.0 

2.6 

18.6 

14.4 

n3 

0.0 

5.3 

16.8 

16.0 

nh2 

0.0 

1.8 

17.1 

13.0 

no2 

0.0 

2.7 

4.9 

R"= 

I-B 

II-B 

III-B 

IV-B 

H 

0.0 

1.9 

19.1 

14.7 

F 

0.0 

5.0 

22.6 

16.8 

CN 

0.0 

1.3 

19.2 

14.6 

n3 

0.0 

2.9 

15.4 

11.0 

nh2 

0.0 

3.5 

16.3 

11.7 

no2 

0.0 

2.5 

17.8 

15.6 

R"'= 

I-C 

II-C 

III-C 

IV-C 

H 

0.0 

1.9 

19.1 

14.7 

F 

0.0 

2.3 

22.1 

14.0 

CN 

0.0 

1.5 

19.9 

14.9 

n3 

0.0 

1.6 

20.1 

9.6 

nh2 

0.0 

1.0 

18.2 

9.4 

no2 

0.0 

1.1 

13.5 

3.3 
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Table  7:  MP2/6-311G(d,p)  average  energy  differences  between  substitution  at  N  and  C  on  the 
tetrazolium  ring. 


Average  Values  (Ec  -  EN): 

F 

-56.9  kcal/mol 

CN 

-20.1  kcal/mol 

n3 

-30.0  kcal/mol 

nh2 

-27.9  kcal/mol 

no2 

+6.0  kcal/mol 

Table  8:  Heats  of  formation  (kcal/mol),  from  G2  theory,  for  parent  isomers. 


Calculated  G3  Heats  of 
Formation3 

Calculated  G2  Heats  of 
Formationb 

I 

248.0 

245.1 

II 

247.1 

246.9 

III 

263.7 

264.3 

IV 

262.8 

259.8 

aG3  Heats  of  Formation  are  taken  from  Satchell  and  Smith13.  bFrom  this  work. 


This  nape  is  Distribution  A:  annroved  for  nublic  release:  distribution  unlimited. 


20 


Table  9:  Changes  in  heats  of  formation  (kcal/mol),  from  G2  theory,  due  to  substituents,  relative 
to  parent  compounds.  Refer  to  Figure  3  for  structures. 


R  — 

I 

r 

average 

R"= 

i  =  r 

R"= 

I 

I’ 

average 

H 

0.0 

0.0 

0.0 

H 

0.0 

H 

0.0 

0.0 

0.0 

F 

17.1 

17.8 

17.5 

F 

-37.5 

F 

19.9 

20.6 

20.2 

CN 

70.7 

71.7 

71.2 

CN 

51.0 

CN 

70.5 

71.5 

71.0 

N3 

108.6 

111.4 

110.0 

N3 

79.7 

N3 

109.9 

112.8 

111.3 

NH2 

18.2 

18.9 

18.6 

NH2 

-10.0 

NH2 

20.7 

21.4 

21.0 

N02 

31.1 

18.1 

24.6 

N02 

18.0 

NO2 

32.8 

19.7 

26.3 

R  — 

II 

IF 

average 

R"= 

II  =  II’ 

R"= 

II 

II’ 

average 

H 

0.0 

0.0 

0.0 

H 

0.0 

H 

0.0 

0.0 

0.0 

F 

20.3 

21.0 

20.6 

F 

-36.0 

F 

20.3 

21.0 

20.7 

CN 

70.2 

71.2 

70.7 

CN 

51.7 

CN 

70.1 

71.1 

70.6 

N3 

109.6 

112.5 

111.0 

N3 

83.1 

N3 

109.6 

112.5 

111.1 

NH2 

19.9 

20.6 

20.3 

NH2 

-10.0 

NH2 

19.8 

20.5 

20.1 

N02 

32.0 

18.9 

25.4 

N02 

18.6 

NO2 

32.0 

18.9 

25.4 

R  — 

III 

nr 

average 

R"= 

hi  =  nr 

R”'= 

III 

III’ 

average 

H 

0.0 

0.0 

0.0 

H 

0.0 

H 

0.0 

0.0 

0.0 

F 

20.5 

21.2 

20.8 

F 

-38.4 

F 

22.8 

23.5 

23.1 

CN 

70.8 

71.8 

71.3 

CN 

50.3 

CN 

71.2 

72.3 

71.8 

N3 

104.8 

107.6 

106.2 

N3 

77.3 

N3 

110.8 

113.7 

112.2 

NH2 

15.4 

16.0 

15.7 

NH2 

-12.1 

NH2 

19.8 

20.4 

20.1 

N02 

N02 

16.6 

NO2 

27.1 

14.0 

20.6 

R  — 

IV 

IV’ 

average 

R"= 

IV  =  IV’ 

R’”= 

IV 

IV’ 

average 

H 

0.0 

0.0 

0.0 

H 

0.0 

H 

0.0 

0.0 

0.0 

F 

19.2 

19.9 

19.5 

F 

-39.6 

F 

19.2 

19.9 

19.6 

CN 

70.7 

71.7 

71.2 

CN 

50.7 

CN 

70.7 

71.7 

71.2 

N3 

104.9 

107.8 

106.3 

N3 

81.0 

N3 

104.8 

107.7 

106.3 

NH2 

15.3 

15.9 

15.6 

NH2 

-11.7 

NH2 

15.4 

16.1 

15.7 

N02 

21.3 

8.3 

14.8 

N02 

18.9 

N02 

21.4 

8.3 

14.9 
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Table  10:  Heats  of  proton  transfer  reaction  (kcal/mol)  for  anion  =  N03‘,  C104' ,  and  N(N02)2’  for 
deprotonation  at  H’,  H”,  &  H”’.  (See  Figure  3) 


Cation 

anion 

H' 

H" 

H’" 

I 

N03- 

-124.8 

-86.1 

-128.5 

C104- 

-102.6 

-63.9 

-106.3 

N(N02)2- 

-112.0 

-73.3 

-115.7 

II 

N03- 

-126.6 

-103.7 

same  as  H' 

C104- 

-104.4 

-81.5 

N(N02)2- 

-113.9 

-90.9 

III 

N03- 

-143.7 

-85.9 

-147.4 

C104- 

-121.5 

-63.7 

-125.3 

N(N02)2- 

-131.0 

-73.1 

-134.7 

IV 

N03- 

same  as  H' 

-68.9 

-142.9 

C104- 

-46.7 

-120.8 

N(N02)2- 

-56.2 

-130.2 
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Figure  1 :  Imidazolium  (a),  triazolium  (b),  and  tetrazolium  (c)  cations. 

Figure  2:  Nitrate,  perchlorate  and  dinitramide  anions  (Oxygen=red,  nitrogen=blue, 
chlorine=green). 

Figure  3:  Resonance  structures  of  isomers  I-IV. 

Figure  4:  The  open  chain  isomer,  azido  formidinium,  V. 

Figure  5:  MP2/6-31  lG(d,p)  IRC  for  proton  transfer  between  isomer  I  and  isomers  II,  III, 
and  IV,  kcal/mol. 

Figure  6:  Isomers  I-IV  with  calculated  bond  orders  in  red. 

Figure  7 :  Isomers  with  a  single  substitution  on  a  nitrogen  or  a  carbon. 

Figure  8:  Isodesmic  reactions  for  all  resonance  structures. 

Figure  9:  Optimized  structures  (A)  and  MP2  and  B3LYP  (in  parentheses)  relative 
energies  (kcal/mol)  for  proton  transferred  structures  that  result  from  optimization  of 
initial  geometries  of  a  tetrazolium  cation  paired  with  a  N03"  anion.  MP2  heats  of  reaction 
(in  brackets)  in  kcal/mol  for  the  reaction:  Tetrazolium  +  N03"  ->  [Tetrazole][H-N03]. 

Figure  10:  Optimized  structures  (A)  and  MP2  and  B3LYP  (in  parentheses)  relative 
energies  (kcal/mol)  for  proton  transferred  structures  that  result  from  optimization  of 
initial  geometries  of  a  tetrazolium  cation  paired  with  a  C104"  anion.  MP2  heats  of 
reaction  (in  brackets)  in  kcal/mol  for  the  reaction:  Tetrazolium  +  C104"  ->  [Tetrazole][H- 
CIOJ. 

Figure  1 1 :  Optimized  structures  (A)  and  MP2  and  B3LYP  (in  parentheses)  relative 
energies  (kcal/mol)  of  the  tetrazolium  cation  paired  with  a  N(N02)2"  anion.  MP2  heats  of 
reaction  (in  brackets)  in  kcal/mol  for  the  reaction:  Tetrazolium  +  N(N02)2"  ->  [Tetrazole] 
[HN(N02)2]  or  (see  text)  [Tetrazolium]  [N(N02)2]. 

Figure  12:  Three  different  geometries  of  dinitramic  acid. 

Figure  13:  The  MP2/6-31+G(d)  potential  energy  surface  for  proton  transfer  from  the 
1,4,5-tetrazolium  cation  to  the  dinitramide  anion  to  form  1,5-tetrazole  and  dinitramic 
acid.  The  relative  energies  including  ZPE  corrections  are  in  parenthesis.  Relative 
energies  in  kcal/mol. 
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Figure  1:  Imidazolium  (a),  1 ,2,4-triazolium  (b),  and  tetrazolium  (c)  cations. 


Figure  2:  Nitrate,  perchlorate  and  dinitramide  anions  (Oxygen=red,  nitrogen=blue, 
chlorine=green). 
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Figure  3:  Resonance  structures  of  isomers  I-IV. 
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•~>v 

Figure  4:  The  open  chain  isomer,  azido  formidinium,  V. 
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Figure  5:  MP2/6-31  lG(d,p)  IRC  for  proton  transfer  between  isomer  I  and  isomers  II,  III,  and  IV, 
kcal/mol. 


This  nape  is  Distribution  A:  annroved  for  nublie  release:  distribution  unlimited. 


26 


This  nape  is  Distribution  A:  unproved  for  nublic  release:  distribution  unlimited. 


27 


IV- A 


) 

-A 


H 


/N\ 

C  N 

W  J/ 

N — N 

\ 


II-A 


H 


/N\ 

H — C  N 

W+  // 

N — N 

/ 

R 


III-A 

H 

H — - R 

w  // 

N — N 


H — 


/Nv 

^N- 

\  / 

N=N 

\ 


H 


) 

[-B 

H 

II-B 

H 

III-B 

H 

IV-B 

— C  N 

W  J/ 

N — N 

\ 

H 

R — C  N 

W+  // 

N — N 

/ 

H 

R — - H 

w  // 

N — N 

/N\ 

R — CT  +N - H 

\  / 

N=N 

\ 

H 

i 

C 

II-C 

III-C 

IV-C 

R 

R 

R 

-C  N 

W  J/ 

N — N 

\ 

H 

/nn 

H — C  N 

W+  // 

N — N 

/ 

H 

H  C^N  H 

\\  // 

N - N 

.N 

H — cr  +  N - R 

\  / 

N=N 

\ 

H 

Figure  7 :  Isomers  with  a  single  substitution  on  a  nitrogen  or  a  carbon. 
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7igure  8:  Isodesmic  reactions  for  all  resonance  structures. 
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Figure  9:  Optimized  structures  (A)  and  MP2  and  B3LYP  (in  parentheses)  relative  energies  (kcal/mol) 
for  proton  transferred  structures  that  result  from  optimization  of  initial  geometries  of  a  tetrazolium 
cation  paired  with  a  NOs"  anion.  MP2  heats  of  reaction  (in  brackets)  in  kcal/mol  for  the  reaction: 
Tetrazolium  +  N03"  ->  [Tetrazole][H-N03]. 
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Figure  10:  Optimized  structures  (A)  and  MP2  and  B3LYP  (in  parentheses)  relative  energies 
(kcal/mol)  for  proton  transferred  structures  that  result  from  optimization  of  initial  geometries  of  a 
tetrazolium  cation  paired  with  a  C104"  anion.  MP2  heats  of  reaction  (in  brackets)  in  kcal/mol  for  the 
reaction:  Tetrazolium  +  C104  ->  [Tetrazole][H-C104], 
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Figure  1 1 :  Optimized  structures  (A)  and  MP2  and  B3LYP  (in  parentheses)  relative  energies  (kcal/mol)  of 
the  tetrazolium  cation  paired  with  a  N(N02)2"  anion.  MP2  heats  of  reaction  (in  brackets)  in  kcal/mol  for 
the  reaction:  Tetrazolium  +  N(N02)2  ->  [Tetrazole]  [HN(N02)2]  or  (see  text)  [Tetrazolium]  [N(N02)2]. 


This  nape  is  Distribution  A:  annroved  for  nublie  release:  distribution  unlimited. 


32 


(+5.2) 

+8.1 

[-117.3] 


(+3.8) 

+6.0 

[-119.3] 


(+8.5) 

+6.7 

[-118.6] 


(+13.1) 

+17.8 

[-106.4] 


(+8.2) 

+7.5 

[-126.4] 


(+4.4) 

+4.5 

[-134.3] 


( - ) 

+0.1 

[-138.6] 


This  nape  is  Distribution  A:  unproved  for  nublic  release:  distribution  unlimited. 


33 


15. 


tV  i* 

2.207  ^ 


(+0.4) 

0.0 

[-138.8] 


(5.3) 

+6.9 

[-135.9] 


(+6.6) 

+8.8 

[-134.0] 


(+9.8) 

+8.5 

[-134.4] 


(+15.1) 

+11.7 

[-131.2] 


This  nape  is  Distribution  A:  unproved  for  nublic  release:  distribution  unlimited. 


34 


Y^» 

Yi* 

YY 

Figure  12:  Three  different  geometries  of  dinitramic  acid. 
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Figure  13:  The  MP2/6-31+G(d)  potential  energy  surface  for  proton  transfer  from  the  1,4,5- 
tetrazolium  cation  to  the  dinitramide  anion  to  form  1,5-tetrazole  and  dinitramic  acid.  The  relative 
energies  including  ZPE  corrections  are  in  parenthesis.  Relative  energies  in  kcal/mol. 
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